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We have observed several cloud cores in the Orion B (L1630) molecular cloud 
in the 2-1 transitions of C^^O, €^"^0 and ^^C^^O. We use these data to show that a 
model where the cores consist of very optically thick C^^O clumps cannot explain their 
relative intensities. There is strong evidence that the C-'^^O is not very optically thick. 
The CO emission is compared to previous observations of dust continuum emission 
to deduce apparent molecular abundances. The abundance values depend somewhat 
on the temperature but relative to 'normal abundance' values, the CO appears to 
be depleted by about a factor of 10 at the core positions. CO condensation on dust 
grains provides a natural explanation for the apparent depletion both through gas- 
phase depletion of CO, and through a possible increase in dust emissivity in the cores. 
The high brightness of HCO+ relative to CO is then naturally accounted for by time- 
dependent interstellar chemistry starting from 'evolved' initial conditions. Theoretical 
work has shown that condensation of H2O, which destroys HCO+, would allow the 
HCO+ abundance to increase while that of CO is falling. 
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1 INTRODUCTION 

There is considerable evidence that protostellar collapse oc- 
curs from an interstellar envelope (scale ~10^ AU) onto, most 
likely, a circumstellar disc (scale ~100 AU) . For a long time 
the collapse is likely to be isothermal and at a low tempera- 
ture (about 10 K). This temperature is lower than the freeze- 
out temperatures of common interstellar molecules (CO 15- 
17 K, NH3 50-60 K, H2O 90 K; see Nakagawa 1980), and, as 
collapse proceeds, the freeze-out timescale becomes shorter 
than the collapse timescale. The expectation is for a long 
isothermal phase at ~10 K when the molecules freeze-out 
onto grains followed by re-heating when they come off again. 
These processes have been modelled by, inter alia, Nejad, 
Williams & Charnley (1990), Rawlings et al. (1992), Bergin 
& Langer (1997) and Charnley (1997). 

There is observational evidence for depletion in cores 
(by about an order of magnitude) in molecules such as CO 
and CS (e.g. Bergin et al. 2001; Tafalla et al. 2002). Most of 
these results have come from observations of nearby molecu- 
lar cloud cores in which low-mass stars are forming in small 
groups (e.g. Jorgensen, Schoier & van Dishoeck 2002). In the 
more distant Orion region where our targets lie, NGC2024 
(in L1630) contains several compact cores which show ev- 



idence for depletion of CO and other molecules (Mauers- 
berger et al. 1992). 

The interpretation as depletion has been criticised, how- 
ever, by Chandler and Carlstrom (1996), who have provided 
evidence that some of the molecular gas is hot. Their crit- 
icisms include: (a) the dust opacity law in the cores will 
be non-standard, (b) dust emission may become optically 
thick at the short submm wavelengths - giving an 'artifi- 
cially low' flux which is interpreted as a low temperature in 
an 'optically thin' interpretation, and (c) the cores may con- 
tain unresolved components which are highly optically thick 
in molecular lines - but thin in submillimetre dust emission. 
Thus the true molecular abundance may be much higher. 
Mangum et al. (1999) also provide evidence that the gas in 
the cores of NGC2024 is hot (>40 K) which would accom- 
modate a normal abundance for CO. 

One problem with NGC2024, however, is that it is near 
strong external heat sources and HIT emission and it is not 
clear that the line emission from the molecular gas arises 
from within the cores from which the dust emission em- 
anates. 

In this paper we report new observations of carbon 
monoxide isotopomers C^^O, C^^O and ^''C^^O from se- 
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Figure 1. LBS17 (22 arcsec beam). Contours are C^^O J=2— >1 
emission integrated from LSR=7 km to 13 km s"'^. Contours 
at 3, 4, 5, 6 and 7 K-km s~^, and peak intensity of 7.60 K-km s~^. 
Greyscale is 850-/im emission with a peaJj flux of 1.76 Jy/beam. 

lected positions near and axound cold cores in the same 
molecular cloud, L1630. These positions have been chosen 
with reference to previously published maps of HCO^, C^*0 
and dust continuum emission from the condensations LBS23 
(also known as HH24-26), LBS17, and LBS18 (e.g. Gibb et 
al. 1995; Gibb & Little 1998, 2000; Phillips, Gibb & Lit- 
tle 2001. These references will be quoted as GLHL, GL98, 
GLOO, and PGL respectively). GL98 and GLOO deduced that 
there was depletion in several cores. The aim of this work 
is to confirm its significance and to evaluate its extent. We 
seek to determine optical depths from the CO isotopomers 
and use the previously published dust emission data together 
with a range of possible emission laws to examine the rela- 
tive abundance. 

In Sections 2 and 3 we justify our need to observe iso- 
topomers as rare as ^^C^**0. In Section 4 we describe the 
observations, which we analyse in Section 5 to deduce the 
depletion in and around the cores. The chemical significance 
of the results is outlined in Section 6. 



2 PROVING DEPLETION? 

Maps of two of the largest cores in the Orion B molecular 
cloud, LBS17 and LBS23, made in high excitation HCO"'" 
and C^*0 (J=2-^l) lines , show striking differences be- 
tween the structures traced by the two molecules. Particu- 
larly apparent is the enhanced brightness of J=3— >2 HCO"*" 
relative to J=2^1 C^^'O in the protostellar core LBS17-H 
(Fig. 2 of GLOO), and of ,7=4^3 HCO+ to J=2^1 C^®0 
in HH25MMS (Fig. 1 of GL98). Figs. 1, 2 and 3 show the 
relation between the C^^O and dust emission. C^^O emis- 
sion is surprisingly weak in some of the clumps which stand 
out prominently in the other species. An interpretation as- 
suming a simple source structure, optically thin C^*0, and 
the dust emissivity law proposed by Hildebrand (1983), im- 
plies widespread reductions of the C^®0 abundance in these 
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Figure 2. HH25MMS (22 arcsec beam). Contours are C^^O 
J=2-^l emission integrated from LSR=7 km s^^ to 13 km s^^. 
Contours at 1.65, 3.3 and 4.95 K-km s^^, and peak intensity of 
5.51 K-km s~^. Greyscale is 850-/um emission with a peak flux of 
1.10 Jy/beam 
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Figure 3. HH24MMS (22 arcsec beam). Contours are C^^O 
J=2— >1 emission integrated from LSR=7 km s"'^ to 13 km s~^. 
Contours at 4.4, 3.2, 4.0, 4.8, 5.6 and 6.4 K-km s~^, and peak 
intensity of 7.66 K-km s~^. Greyscale is 850- /im emission with a 
peak flux of 2.21 Jy/beam 
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clumps by factors of 13 to 57 compared to what we take 
to be its canonical value of 2xl0~^ (GL98 and Frerking, 
Langer & Wilson 1982). Most of the clumps appear to be 
bound objects of several solar masses (for LBS23 at least). 

GL98 consider other ideas to explain the apparent re- 
duction in abundance for the clumps in LBS23. Abnormal 
dust properties do not appear to be required because dust- 
emission- derived masses for the clumps are in good agree- 
ment with their virial masses. 

Another possibility is the absorption of C^**0 emis- 
sion from intervening low excitation temperature gas. If this 
were the case then the effect would be seen as strong self- 
absorption in the optically thicker CO lines. However, there 
is no evidence for such self-absorption in the CO spectra 
(GH93). Since the observed antenna temperatures measure 
the excitation temperature where the optical depth becomes 
equal to unity, a comparison of the CO and C^^O antenna 
temperatures would then suggest that the cloud cores are 
externally heated (see e.g. GL98). 

IftheC^**0 emission arises in optically thick, unresolved 
sub-clumps the need for depletion might be avoided as all 
the C^^O might not be detected. Accordingly we have made 
careful observations of weak isotopomers, including C^**0, 
C^'^O and ^^C^®0, to seek to eliminate a sub-clump model. 
To illustrate how this can be achieved we consider a simple 
two-component model in the next section. 



3 OPTICALLY THICK CLUMPING MODEL 

In Section 5 of GL98 the authors describe a two-component 
model ('envelope' and 'sub-clump') where the envelope 
(component 1) is optically thin (depth ri, linewidth Aui) 
in the less abundant CO isotopomers, such as C^*0, and 
the sub-clump material (component 2) is optically thick 
(depth T2, linewidth /S.V2) in these same lines. Extending 
their analysis it is also possible to define an apparent deple- 
tion factor which describes how the presence of the optically 
thick clump mimics depletion within the beam in the sim- 
plistic interpretation as a single component. The CO abun- 
dance for such a source would be proportional to (TXiAui 
-\- TX2Aii2r2). However, in the standard optically thin anal- 
ysis, not knowing T2, we would assume the CO abundance 
is proportional to (TJi Aui -I- TX2^V2). We will thus deduce 
CO abundances that are too low by a ratio 

niAvi+Tl2Av2 _ R12 + I 
TJiAvi -FTXjAiJara R12 + T2 ^' 

where R12 is defined in GL98 as the ratio of the integrated 
intensities of components 1 and 2. The optical depth of the 
sub-clump component T2 in the denominator causes A to 
be less than unity, and demonstrates how an optically thick 
sub-clump can mimic depletion. If we assume that the abun- 
dance is in fact normal then we can use our model along with 
estimates of R12 and A to deduce information concerning the 
nature of the sub-clumps. 

For example, the clump HH25MMS (GLHL, GL98 and 
PGL) stands out very prominently in dust continuum and 
4^3 HCO^ emission but it is barely detectable as an in- 
dependent source on the north-south ridge of 2^1 C^^O 
emission. The integrated intensity on the ridge is typically 
2.5 Kkms-i while HH25MMS adds about 0.8 Kkms'^ to 



it. Thus R12 = 2.5/0.8 — 3.1. Also the apparent depletion 
factor, A, is about 1/30 (according to GL98). Hence, from 
equation Q T2 — 120. 

GL98 showed that the beam filling factor of the sub- 
clumps was given by 



i?i2Aw2 -I- n Awi 

The 4^3 HCO"*" emission is closely correlated with the 
continuum dust emission, rather than the C^*0, but has 
a linewidth similar to the C^*0 so Avi = Av2. Also for 
HH25MMS ri ~ 0.2 so that the sub-clump filling factor is 
fc ~ 0.06. These would be sub-clumps with a very high C^^O 
optical depth and a very small beam filling factor in the 22 
arcsec JCMT beam. 

We shall assume that the JCMT beam, of radius R, con- 
tains n identical sub-clumps of radius r. Then fc — n{r/R)'^. 
We know n > 1 since HH25MMS is resolved. It follows that 
r < 2.4 arcsec or 1000 AU at the distance of L1630 (400 
pc). Since the apparent depletion is high, most of the mass 
within the beam at HH25MMS is in the sub-clumps. GLHL 
quote a mass of 8 M© in the core of HH25MMS. The sub- 
clump mass is then ~ (8/n) Mq and its density is greater 
than 2.5x10* cm~^. The density of the envelope material 
is much lower. Its mass is l/30th that of the clumps, and 
taking the mass to be within a beamwidth gives a density 
5x 10'* cm~^. In summary, if we seek to retain normal abun- 
dances by assuming sub-clumping it would appear that the 
sub-clumps must have radii less than 1000 AU and densities 
greater than 2.5x10® cm~^. 

As described by GL98, the most obvious problem with 
the 'sub-clump' interpretation seems to lie in the fact that 
the HCO+ intensity follows the submillimetre dust contin- 
uum emission well (rather than the envelope observable in 
C'^^O) yet its brightness temperature is so high that it must 
have a beam filling factor near unity, which is much greater 
than that deduced for the sub-clumps from C**0. Nonethe- 
less it is important to test the structure within the clumps 
either via interferometric observations or beam-matched ob- 
servations of weaker CO isotopomers. 

We may predict C^'^O line intensities, assuming r2=120. 
The component 1 intensity should decrease by a factor of 
5.4 compared to C**0, since it is optically thin, but that 
from component 2 will be unaltered as its optical depth will 
stiU be extremely high (=120/5.4). Referring to Tabled] the 
intensity will be 

7a = 7ai +TI2 = 0.34 K + 0.2 K = 0.54 K (3) 

so the ratio of the C^^O to C^'^O intensities wiU be 2.6. This 
might easily be misinterpreted as indicating that the C'^^O 
was just slightly optically thick. Observations of even rarer 
isotopomers - so that the optically thin 'envelope' contribu- 
tion becomes negligible while the optically thick core com- 
ponent still remains - are required. Table Ogives predicted 
line intensities for the 'depletion' and 'sub-clump' models. It 
can be seen that observations of '^^C'^^O are clearly capable 
of distinguishing between them. 
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Table 1. Line optical depth and brightness (K) predictions for the 'sub-clump' and 'depletion' models, and observed values, for 
HH25MMS. 



Molecule 


Sub-clump 
optical depth 


Sub-clump model 

(K) (K) 

sub-clump Envelope 


Total 
(K) 


Depletion 
model 
(K) 


Observed 
(K) 




120 


0.34 


1.1 


1.4 


1.4 


1.3 


c"o 


22 


0.34 


0.2 


0.55 


0.25 


0.28 


13C180 


2 


0.34 


0.02 


0.35 


0.025 


<0.08 



4 OBSERVATIONS 
4.1 CO observations 

Observations of the J=2^1 transitions of C^*0, and 
isqISq .^^gj.g carried out over the period of 1998 February 
16-28 at the James Clerk Maxwell Telescope, using the com- 
mon user SIS receiver A2 (Davies et al. 1992) and the digital 
autocorrelation spectrometer DAS. The original mixer in A2 
had been replaced by one from NRAO. The C^'^O observa- 
tions were made assuming a rest frequency of 224.714 GHz, 
the C^*0 observations at 219.560 GHz and the "G^**0 ob- 
servations 209.419 GHz. Key positions were selected from 
our previous maps of the cores, and, in addition, a new map 
of C^**0 in LBS18 was made. The positions chosen are listed 
in Table 121 The positions of cores were observed, as well as 
a selection of references off core where there was sign of 
continuum emission. 

The pointing and antenna temperature were measured 
with reference to the nearby source OMCl. The pointing was 
consistent typically to within about 3 arcsec. On 16th and 
17th February the apparent antenna temperatures of lines 
from the frequently observed calibrating source OMCl var- 
ied with time when observing the C^^O and ^'^C^^O lines. 
The intensity of the C^^O line varied systematically with 
time between 0.51 and 0.80 of that of the Observatory ref- 
erence spectrum for OMCl C^^O, Tl = 7.5 K. The C^^O 
intensities we quote have been systematically corrected for 
this variation, which brought three of the observations into 
agreement with similar ones described by GL98. 

There was no appropriate Observatory reference for 
that part of the passband containing ^'^C^^O or C^^O. Al- 
though ^'^C^^O was not clearly discernible for OMCl on 16th 
February the intensities of several bright lines in the pass- 
band containing it varied systematically with time between 
0.42 and 0.64 of their value on 18th and 19th February. It 
has been assumed that the values for 18th and 19th Febru- 
ary are correct, and our ^^C^**0 intensities scaled on this 
assumption. The C^^O line observed in OMCl showed no 
significant variation between 17th, 18th and 20th February, 
yielding TX ~ 1.9±0.1 K. This value has been assumed cor- 
rect, and no scaling applied to our observed C^^O intensities. 

An observation of CO in OMCl on 17th February 
yielded a line of brightness 0.9 of the observatory refer- 
ence value. The reason for the time variation of C^*0 and 
isqIBq jjjjg intensities on 16th and 17th February has not 
been established. Based on the authors' experience, possibil- 
ities might include incorrect adjustment of one of the tele- 
scope mirrors, problems in SIS receiver calibration (see, e.g. 
Davies et al. 1992) or uncertainty in receiver sideband cal- 



Table 2. Positions studied in this investigation (B1950). 



Source 


h 


RA 
m 


s 


O 


DEC 


ff 


Core? 


HH24MMS 


05 


43 


34.7 


-00 


11 


49.0 


yes 


HH25MMS 


05 


43 


33.8 


-00 


14 


45.0 


yes 


LBS17H 


05 


43 


57.1 


-00 


03 


44.3 


yes 


LBS18S 


05 


43 


54.2 


+00 


18 


23.0 


yes 


LBS23E1 


05 


43 


33.8 


-00 


14 


36.0 


no 


LBS17F 


05 


43 


54.2 


-00 


02 


48.0 


no 


LBS18A 


05 


43 


55.2 


-fOO 


18 


57.0 


no 


LBS18B 


05 


43 


54.7 


+00 


18 


48.0 


no 


LBS23A 


05 


43 


35.0 


-00 


11 


00.0 


no 



ibration, though the latter might not be expected to vary 
with time. The receiver was tuned for double sideband per- 
formance and it is expected that the sideband ratio should 
be no worse than 1.2. Changing the LO frequency to observe 
the same C^*0 line in both sidebands on the 17th, when the 
observed strength was low (though not a measurement of the 
sideband ratio) produced very similar line intensities, tend- 
ing to suggest that the sideband ratio was not responsible. 

Fig. |1] shows the new map of the LBS 18 region along 
with the 850 map from PGL, and Figs. |3 to |H| show the 
new spectra obtained for HH24MMS, HH25MMS, LBS17H 
and LBS18S. 

4.2 Continuum observations 

The continuum observations used in this paper were made 
using the submillimetre bolometer array receiver SCUBA 
and are described in PGL. For the purposes of this paper 
we have smoothed their 850/im maps to the same resolution 
as that of our CO data, i.e. 22 arcsec. 



5 RESULTS AND ANALYSIS 

Table 13 contains the data taken from the new observations. 
The antenna temperatures, TX, derived from the CO ob- 
servations for all three species at each of the positions are 
shown. The observations were made with a channel spac- 
ing and noise bandwidth of 0.45 kms~^. Four channels 
were averaged to give the temperatures quoted in Table 
13 (Note that the spectra are show in the original spec- 
tral resolution of 0.45 kms~^.) The integrated intensity 
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Figure 4. LBS18 (22 arcsec beam). Contours are C^^O J=2^1 
emission integrated from LSR=7 km s"'^ to 13 km s~^. Contours 
at 1.25, 1.75, 2.25, 2.75 and 3.25 K-km s"!, and peak intensity 
of 3.35 K-km s~^. Greyscale is 850-/im emission with a peak flux 
of 0.65 Jy/beam. 
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Figure 5. HH24MMS: C^^^O, C^^O and ^^C^^O spectra at the 
core position (RA=05''43'"34.7* DEC=-00° 11'49.0"). 
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Figure 6. HH25MMS: C^^^O, C^^'O and "C^^O spectra at the 
core position (RA=05''43'"33.8^ DEC=-00° 14'45.0"). 



7int(Kkms~^), which is derived from the antenna tempera- 
ture (/int=/ T^dv, integrated over ±1.79 kms~^), for each 
species is also shown. 

The Astronomical Image Processing Software (aips) 
package was used to analyse the continuum data. The 850 
^m and 450 /xm continuum emission maps were convolved 
to an effective resolution of 22 arcsec. The fluxes F (Jy) for 
each position listed in Table|21were then taken from the con- 
volved maps and are listed in Table O The uncertainty in 
these values is 0.02-0.03 Jy. 
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Table 3. C^^O, C^^O and ^'^C^^O line temperatures and integrated intensities, and 850 fira flux measurements at each position. Note 
that the noise levels quoted are derived from spectra binned by four channels and are thus a factor of 2 lower than might be expected 
from a direct comparison with Figs.|3to|Sl 





C 


18o 


C 




13q18 


O 


850 /im 


Source 




hnt 


Ta* 


/int 


Ta' 


lint 


F 




(K) 


(K km s-l) 


(K) 


(K km s-l) 


(K) 


(K km s-i) 


(Jy) 


HH24MMS 


2.8±0.03 


5.8±0.07 


0.67±0.01 


1.75±0.03 


0.057±0.020 


0.09±0.03 


4.76 


HH25MMS 


1.3±0.03 


2.5±0.09 


0.28±0.01 


0.80±0.02 


<0.080±0.045 


<0.02±0.05 


2.36 


LBS17H 


2.6±0.03 


6.0±0.07 


0.64±0.01 


1.80±0.03 


0.088±0.048 


0.16±0.06 


3.79 


LBS18S 


1.5±0.03 


3.0±0.16 


0.37±0.01 


0.80±0.03 


<0.072±0.048 


<0.02±0.06 


1.34 


LBS23E1 


1.8±0.06 


3.4±0.15 


0.30±0.02 


0.80±0.05 






1.97 


LBS17F 


3.9±0.06 


7.3±0.15 


0.90±0.02 


2.20±0.05 






1.85 


LBS18A 


2.6±0.06 


3.0±0.15 


0.45±0.02 


0.75±0.05 






0.82 


LBS18B 


3.0±0.06 


3.4±0.13 


0.42±0.02 


0.90±0.05 






0.85 


LBS23A 


2.6±0.10 


6.0±0.23 






<0.05±0.015 


<0.12±0.02 


1.08 
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Figure 7. LBS17H: C^^^O, C^^O and "C^^^O spectra at the core 
position (RA=05''43™57.1'' DEC=-00°03'44.3"). 



1 ' 








13C180 


, 1 , 



-20 20 

Velocity / (km/s) LSR trame +10.00 km/s (Radio Def'n) 

Figure 8. LBS18S: C^^^O, C^'^O and i^ci^O spectra at the core 
position (RA=05''43™54.2'' DEC=+00° 18'23.0"). 



5.1 Temperatures and masses 

From the fluxes at 850 /im and 450 we seek to derive the 
temperatures and masses. The results depend on the opacity 
law, which is uncertain. Rather than seeking to justify a law 
we chose to investigate the effect of varying it by compar- 
ing the different results obtained assuming two typical but 
reasonably extreme examples: the laws proposed by Hilde- 
brand (1983) (/3=2) and by Testi & Sargent (1998) (/?=1.1). 
These predict very similar opacities at 850 /xm, but different 
temperatures from 850/450 flux ratios. We accordingly de- 
rive temperatures from both these laws, and compare them 
with already published temperatures from grey body spec- 
tral flts and ammonia inversion line observations to deduce 
plausible temperature limits. These are then used with the 
850 fxm opacities to deduce masses. 

If the absorption coefficient 1^(1/) is written in terms of 



the dust and molecular hydrogen density as 
— KdPd = «:H2''^H2 



(4) 



the total number of hydrogen molecules N{H2) implied by 
observation of a flux F^, from dust assumed to be optically 
thin, is 



iV(H2) = 



(5) 



where D is the distance from Earth. 

The temperature of the gas can be derived by forming 
the ratio of the submillimetre continuum fluxes at the two 
wavelengths and comparing it with that predicted for dif- 
ferent temperatures using equation |S| Flux-ratio maps were 
produced by using the comb task in AiPS to divide the 850 
/xm (15 arcsec resolution) maps by the 450 /xm (15 arcsec res- 
olution) maps. The ratio, r(850/450), was measured at each 
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Figure 9. Flux-ratio map (22 arcsec) of the LBS18 region plotted 
over the 850-/xm (22 arcsec) continuum map. Contours levels are 
at 0.3, 0.4, 0.5, 0.6, 0.7 and 0.8 ratio, with a peak of 0.9 ratio at 
the dust peak. 



point from these maps and used to find the temperature. The 
flux-ratio map of LBS18, which has been smoothed to a res- 
olution of 22 arcsec, has been plotted over the 850 /im (22 
arcsec) continuum map in Fig. |51 where it can be seen that 
the contours representing the higher values of r (850/450) 
are coincident with the dust peaks. This may be due to (a) 
a lower temperature, (b) a flatter emissivity law, or (c) the 
450 /im emission being optically thick. 

The resulting value of Td for each position can be found 
in Table 2] for both values of 13. As expected, the tempera- 
tures are higher for the lower value of /3. The errors are de- 
rived using the calibration uncertainties and random noise 
errors quoted in PGL. 

Assuming the dust is optically thin can produce mis- 
leadingly low temperatures and high derived masses if the 
dust is in fact optically thick at 450 /im. However PGL used 
a radiative transfer program to model the cores that was not 
limited by this assumption, but derived very similar temper- 



Table 4. Temperatures (K) derived from published grey body 
fits, our observations assuming f} = 1.1 and 2, and NH3 observa- 
tions of Harju et al. 



Source Tgb Ti.i T2 TNH3 



HH24MMS 


20±5 


10 


+4 
-3 


7±1 


15±3 


HH25MMS 


17-35 


38: 


f? 

-20 


13+3 


15±2 


LBS17H 


<20 


14 


+9 
-5 


9+3 
^-2 




LBS18S 




14: 


H2 

-5 






LBS23E1 




37: 


f? 

-22 


13ti 




LBS17F 




16 


+'! 
-8 






LBS18A 




23: 


f? 
-14 






LBS18B 




24: 


-14 






LBS23A 




37: 


f? 
-23 


13t» 





atures to those found here. In reality the results we find rep- 
resent an average over the beam. It is also quite likely that 
there are temperature gradients within the beam, either de- 
creasing towards the centre of the core (e.g., an externally 
heated protostellar cloud) or increasing (e.g., an interstel- 
lar cloud heated by a central young star). If there are, then 
the derived beam-averaged temperature will be weighted 
towards the hotter dust which radiates more strongly. This 
means that the derived dust masses will be too low (colder 
dust being underweighted) . CO column densities are less 
sensitive to temperature, so the overall effect is to derive 
CO abundances which are higher than the true ones - i.e. 
to make the depletion appear less marked than it really is. 

Further information is available which can be used to 
constrain the temperatures. Grey body fits to the spectra of 
dust continuum emission give rd = 17-35 K for HH25MMS 
(Gibb and Davis, 1998), rd=20±5 K for HH24MMS (Ward- 
Thompson et al., 1995) and Td <20 K for LBS17H (GLOO). 
It is also of interest to compare the derived temperatures 
with those obtained from NH3 inversion lines. Harju et al. 
(1993) used NH3(1,1) and (2,2) observations of 40 arcsec res- 
olution to derive the kinetic temperature of 43 star forming 
cores, 16 of which are associated with the Orion L1630 and 
L1641 clouds. They mapped LBS23, finding temperatures of 
about 15K at the positions of HH24MMS and HH25MMS. 
These results and those from the grey body fits are included 
in Table |1| In general the kinetic temperature shows little 
variation over their maps. Their results show that the av- 
erage kinetic temperature is ~15.7 K for the Orion cores. 
Only two of their cores showed temperatures greater than 
20 K. Matthews and Little (1983) also mapped LBS23 with 
130 arcsec resolution; at HH25 they found a temperature of 
14 K which is similar to that of Harju et al. As the cloud 
is extended and their beam contained a substantial region 
surrounding the core it is likely that a temperature of 14 K 
is representative of the envelope surrounding the core. 

On the basis of Table |1| and the general considerations 
above we regard conservative limits for the dust temper- 
atures to be 20±5 K for HH24MMS and HH25MMS and 
15±5 K for LBS17H and LBS18S, while for the off-core 
envelope positions we take T=15±5 K. We then use these 
temperatures and the 850 /xm fluxes together with the 850 
^m opacity of Testi and Sargent to deduce gas masses via 
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equation 5. These are used in Section 5.3 to derive relative 
CO abundances. 



5.2 Carbon monoxide optical depths 

We can see immediately that the observed isotope line ratios 
listed in Table for HH25MMS are in much better agree- 
ment with the predictions of the depletion model than with 
those for optically thick sub-clumps. It should be noted that 
although optically thick undepleted sub-clumps are unlikely 
to explain the apparent depletion the existence of depleted 
clumps is not ruled out. 

Table |3 shows the line temperatures and their ratios 
(X(18/17), denoting the ratio of Taj^ and TX-^.,, the line 
temperatures of C^^O and C^'^O respectively). The mean 
value of X(18/17) in the cores is 4.28; towards the other 
positions 6.1. The latter value is in good agreement with the 
solar expectation of 5.4, suggesting that these positions are 
optically thin, and that the [C^^O]/[C^*0] ratio is normal. 
The former suggests that C^^^O may be slightly optically 
thick in the cores, but certainly not sufficiently to discredit a 
former conclusion that CO is depleted in several of the cores 
of LBS23 (GL98). The [C^'' O\/[^^C^^0] ratios are consistent 
with the solar value. 

The above constitutes good evidence that C^^O and 
^^C^**0 are optically thin. 

The optical depths of the three species in the individual 
positions have been calculated. Firstly, the optical depth 
of C'^^'O can be calculated using equation |S| where Tai7 
is the line temperature of C^'^O and T^-^g is that of C^*0 
and Ti7 and rig are the optical depths of C^'^O and C^^O 
respectively. 



Table 5. Optical depth of C^^O at each position. 



Source 


(K) 


Tt 
(K) 


Jf(18/17) 


Tl -7 
' i / 




HH24MMS 


0.67 


2.8 


4.2 


0.12 


0.02 


HH25MMS 


0.28 


1.3 


4.6 


0.08 


0.04 


LBS17H 


0.64 


2.6 


4.1 


0.13 


0.02 


LBS18S 


0.37 


1.5 


4.1 


0.13 


0.04 


LBS23E1 


0.30 


1.8 


6.0 


-0.05 


0.02 


LBS17F 


0.90 


3.9 


4.3 


0.11 


0.03 


LBS18A 


0.45 


2.6 


5.8 


-0.03 


0.05 


LBS18B 


0.42 


3.0 


7.1 


-0.12 


0.05 


LBS23A 




2.6 


7.3 


-0.13 





Table 6. Temperatures and CO abundances 



Source 



T 
(K) 



(10-8) 



Xci''0 X13C180 
(10-«) (10-10) 



HH24MMS 
HH25MMS 
LBS17H 
LBS18S 

LBS23E1 

LBS17F 

LBS18A 

LBS18B 

LBS23A 



20±5 
20±5 
15±5 
15±5 

15±5 
15±5 
15±5 
15±5 
15±5 



2.4 
1.9 
2.0 
2.9 

1.7 
4.8 
3.6 
3.9 
5.4 



1.1 

0.8 
+0.9 
0.7 
+0.8 
0.7 
1.1 
0.9 

+0.9 
0.8 
+2.6 
2.2 
1.9 
1.6 
+2.1 
1.8 
+2.9 
2.4 



0.6 
0.5 
0.5 
0.6 

0.4 
1.2 
0.9 
1.0 



+0.3 
0.2 
+0.2 
0.2 
+0.3 
0.2 
+0.3 
0.3 



+0.2 
0.2 
+0.6 
0.5 
+0.5 
0.4 
+0.6 
0.5 



2.8 

<1.3 
4.1 
<1.5 



1.6 
1.4 
+2.0 
2.0 
+2.3 
-2.0 
+3.0 
3.0 



— <10. 



D+6.9 
'-6.1 



-5.4T17 



(6) 



By assuming that solar values hold, ris = 5.4ri7, equa- 
tionl^was solved to obtain a value for ri7 at each position. 
The values of ri7 are displayed in Table |^ with the errors 
calculated for each value. Multiplying ri7 by 5.4 gives rig 
at each position. One can also calculate ri3^ig, the optical 
depth of "C^*0 by dividing ri7 by 16.6. It can be seen from 
these results that all three species of CO are optically thin 
in the regions where the abundances have been calculated 
and can therefore be used effectively in the calculations. The 
negative optical depth values in Table |K| are likely to be due 
to noise errors, rather than weak maser emission. 



5.3 Relative abundance of CO 

For an optically thin transition the total number of 
molecules within a beam 6 (FWHM arcsec) from a source 
at D (kpc) is given by 



NiMOV) = 2.5 X 10* 



D_ 

04 



22 



55.5 
B 



Tndv 



(2J + 1) /or" ^ 

J^Fj{l-R) \ M 



(7) 



where B is the rotational constant (GHz), /i the dipole mo- 
ment (debyes) , J the upper level of the transition, Fj is the 



fractional population of level J, J T^dv is the integrated line 
intensity (Kkms-'^), and R is given by 



R 



exp(2|f^) 
exp(- 



(8) 



fcTB 



where Tb is the cosmic microwave background temperature 
(2.7 K). 

The relative abundance of the CO isotopomer can be 
estimated by computing N{MOL) from equation and di- 
viding by A'^(Il2) from equation|^ (To calculate the relative 
abundance, the integrated intensities in Tablej^were divided 
by r?B(= 0.8) to correct for beam efficiency.) 

The C^^O abundance, Xci»0) "^^-^ calculated in the 
same way as for C"0 and "C^^^O but the answers obtained 
for HH24MMS, HH25MMS, LBS17H, LBS18S, and LBS17F, 
were muhiplied by (5.4/X), where X=X(18/17) (see Table 
[SJ, to allow for the finite optical depth of C^^O. The result- 
ing relative abundances of C^'^O, C^^^O and ^^C^^O at each 
position can be found in Table |S| 



6 DISCUSSION 
6.1 Morphology 

It appears, from comparing SCUBA and line results, that 
the cores tend to be found on the edges or ends of fila- 
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ments. There are two kinds of core structure: ones where 
the HCO^ and dust are coincident but C^^O is very faint or 
absent (LBS17H, HH24MMS), others where HCO+, C^^O 
and dust peak successively towards the cloud edge (LBS18S 
and the Launhardt et al. (1996) core in LBS17 (see GLOO)). 
The problem is to distinguish the effects of chemistry and 
excitation. 



6.2 Depletion 

The canonical values for the abundance of C^**0, C^'^O and 
^^C^*0 may be taken to be 2 x 10"'^, 4.7 x 10"* and, assum- 
ing simple proportion, 2.7 x 10~^ respectively, from Frerk- 
ing et al. (1982). Dividing these values by the corresponding 
abundances in each position, one can estimate the deple- 
tion factor, d, for each isotopomer. Table [7| displays these 
depletion factors and one can see that in all cases the CO is 
depleted significantly. 

Relative to the canonical values, the isotopic species 
appear to be depleted by factors of about 10 on the cores, 
and approximately half that in the envelopes. The values 
are typically less than those quoted by GL98, the differences 
being due to the assumption of slightly higher temperatures, 
careful matching of effective beam sizes and correction for 
the optical depths of C^**0. 

A natural explanation for the depletion observed in 
these cores is the condensation of CO molecules on dust 
grams smce the density (-10*^ cm"^: GLHL, GL98) is high 
enough to reduce the freeze-out timescale to less than the 
dynamical timescale. In addition, grain growth by coagula- 
tion is expected in such regions which gives rise to increased 
dust emissivity (Ossenkopf 1993). The consequence of these 
two processes is that dense collapsing cores become very well 
defined in submillimetre continuum emission while simul- 
taneously becoming very poorly defined in C^*0 emission, 
such as we observe here. 

A number of studies of other star-forming regions have 
also arrived at the conclusion that the CO abundance is 
lower than the canonical value. Studies by Kramer et al. 
(1999), Willacy, Langer & Velusamy (1998) and Caselh et 
al. (1999) have found that the C^^O abundance is reduced 
typically by factors ranging from a few to an order of mag- 
nitude towards the lines of sight with highest H2 column 
density. The study of Kramer et al. (1999) shows that the 
highest reduction in the C^*0 abundance occurs in the cold- 
est regions. In all cases, these authors also conclude that 
freeze out of CO onto cold dust grains is responsible for the 
observed levels of depletion. 

Although this investigation has shown that CO appears 
to be depleted in the cores, one cannot be sure that the fac- 
tors are correct, as we are calculating the abundance using 
averages along the line of sight. If, for example, the cores are 
much colder than the envelopes, there will be a higher de- 
pletion factor in the cores and a lower one in the envelopes. 
This is one area which needs to be studied further in or- 
der to find a way of deriving the emission produced at the 
core positions rather than averaged along the line of sight. 
One must also keep in mind that the canonical values of the 
abundance have been calculated using line of sight observa- 
tions which may contain dust and gas that is a lot less dense 
than that found in protostellar cores so comparing the core 



Table 7. Depletion factors 



Source 


T (K) 






rfl3clSo 


HH24MMS 
HH25MMS 
LBS17H 
LBS18S 


20±5 
20±5 
15±5 
15±5 


c+4.6 
».0_2.7 


c:+4.6 
"'■0-2.6 

9.2+5-0 
10 2+** '5 


21. 3+?'^ 

i8.etf° 


LBS23E1 

LBS17F 

LBS18A 

LBS18B 

LBS23A 


15±5 
15±5 
15±5 
15±5 
15±5 


^•^-1.4 

0.0_2 

o 7+3.1 
■^•'-1.3 


4 0+3-5 
^■'-'-1.4 

0+4.5 
^■•^-1.9 
4(5+3.8 


9 c:+3.2 
^•0-1. 



abundances with these values may not give a true depletion 
factor. 

6.3 Implications for chemical models 

Rawlings et al. (1992) show that as a protostar collapses it 
is plausible that the chemistry occurring while molecules de- 
plete out onto dust grains leads to an increasing HCO+/CO 
ratio (by up to 50), although their absolute abundances fi- 
nally fall. A requirement for this to happen is a high initial 
H2O abundance such as might result from initial conditions 
following a shock through already processed material, but 
not from initial conditions which are of simple atomic form. 

Observations of water have been made towards the re- 
gions encompassing HH24MMS and HH[25MMS with ISO, 
but these only trace a hot dense component in shocked gas 
(Benedettini et al. 2000). Perhaps more reliable but still 
less than ideal are the observations of NGC2024 with the 
Submillimeter-Wave Astronomy Satellite (SWAS) by Snell 
et al. (2000). Snell et al. estimate that the abundance of 
water (strictly only the ortho species) to be 6 x lO"^'', con- 
siderably lower than expected but in agreement with the 
notion of ongoing and widespread depletion by freeze out 
onto dust grains. 

Other models have been developed which have differ- 
ent assumptions. Bergin & Langer (1997) followed chemistry 
during the collapse using a gas-grain chemical code, includ- 
ing desorption and depletion, but starting from an initial 
pristine state with elements only. The increase in HCO+/CO 
suggested by observations was not produced by the model. 
Neither was it produced by the different type of model of 
Charnley (1997). 

Bergin et al. (2000) have proposed an alternative 
method of reducing the CO abundance, prompted by re- 
sults from SWAS. In their model CO abundance can drop 
(even in warm clouds, say, 30 K) as a result of destruction 
by He+ with the oxygen atoms freezing out and becoming 
locked in water ice. Unfortunately their discussion excludes 
HCO+, and it should also be noted that the CO destruction 
occurs on extremely long timescales (greater than a million 
years), much longer than the collapse timescales for these 
cloud cores. 

The model of Rawlings et al. (1992) predicts an abun- 
dance for IICO+ in good agreement with a value we have 
derived for HH25MMS (5 x 10"^°, GLHL). In addition, the 
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modelling of Caselli et al. (2002) shows that the HCO+ 
abundance may be used to estimate the maximum degree 
of depletion. Their Model 3 (shown in their Fig. 8) also pre- 
dicts a value for the HCO"'" abundance close to our LVG 
estimate for HH25MMS. At this value (which is also at a 
radius close to the beam radius for our JCMT observations) 
the depletion factor is ~20, in good agreement with the lev- 
els of CO depletion we derived above. 

However,it should also be noted that from our LVG 
modelling (GLHL) we find that the HCO^ is optically thick, 
and thus any decrease in the HCO"*" abundance in the centre 
of these cores (at radii less than ~3000-4000 AU) is masked. 
Observations and modelling of rarer species will help to an- 
swer this question. We are currently analysing observations 
of HCO"*" and DCO"*" isotopomers on the cores to study this 
phenomenon more carefully. 

Since star-formation has already taken place in L1630 
(e.g. Lada ct al. 1991) the initial conditions for the chem- 
istry will not be pristine. This supports our consideration 
of the model of Rawlings et al. (1992), which starts from 
'evolved' conditions (including a high abundance of water). 
Only this model predicts an actual increase of HCO'*' while 
CO is declining. It thus provides a natural explanation for 
the brightness of high excitation HCO"'' emission relative to 
CO emission on cores. It thus appears that the role of the 
initial conditions is very important in determining the later 
behaviour of the chemistry. 



7 CONCLUSIONS 

We have used new observations of C^^O, C^**0 and ^^C^^O 
along with existing 450 jj,m and 850 ^m continuum observa- 
tions to study four regions in the Orion B molecular cloud. 
This has allowed the determination of optical depths for the 
isotopic species and the elimination of a model by which the 
cores contain many optically thick sub-clumps to explain the 
low isotopic CO emission. 

If the dust emissivity falls within what are presently 
considered to be reasonable limits, CO is depleted in both 
the cores and the surrounding regions, but by a higher fac- 
tor in the cores. Relative to the canonical abundances of 
Prerking et al. (1982), the CO appears to be depleted by 
about a factor of 10 at the core positions and a factor of 
5 in the envelopes. Although the dust observations do not 
allow the temperatures to be constrained for /3=1.1, we be- 
lieve it unlikely, from the NH3 observations of Harju et al. 
(1993) and Matthews and Little (1983), that temperatures 
will be greater than 15-20 K. It seems unlikely that the 
temperatures exceed the dust derived values in the sources. 

The low abundance of CO and the brightness of HCO'*' 
emission in the cores tends to support the model of Rawlings 
et al. (1992) which suggests that the amount of HCO'*' in 
the core increases as that of CO decreases. The HCO^ abun- 
dance and observed levels of CO depletion are also in good 
agreement with a recent model of molecular ion chemistry 
in collapsing cores by Caselli et al. (2002). 
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